WORLD TNTEU-ECrUAL PROPERTY ORGANIZATION 
[ntcmaiiunal Bureau 




PCX 

INTERNATIONAL APPLICATION PUBUSHED UNDER THE PATENT COOPERATION TREATY (PCT) 




(51) International Patent Classification 6 ; 




(11) TnternaUonaJ Publication Number: 


wo 98/20361 


GOIS 7/52 


Al 


(43) Internationa] Publication Date: 


14 May 1998 (14.05.98) 



(21) International Application Number: PCT/US97/20162 

(22) International Filing Date: 4 November 1997 (04.1 1.97) 



(30) Priority Data: 
08i746.360 



8 Noverabcr 1996 (08.11.96) US 



(71) Applicant: RESEARCH CORPORATION TECHNOLOGIES, 

INC. [US/US]; Suite 600, 101 N. WiJmot Road. Tucson, AZ 
85711-3335 (US).' 

(72) Inventor: CHRISTOPHER, Ted; Apartment I, 51 Lilac Drive, 

Rochester, NY 14620 (US). 

(74) Agents: DiGIGLIO, Fi^aiik, S. ct al.; Scully, Scott, Muiphy & 
Pi^ser, 400 Garden City Plaza. Garden City, NY 11530 
(US). 



(81) Designated States: AU, CA, IP, European patent (Al\ BE, 
CH, DE, DK, ES, FI, PR, GB, GR, IE, IT, LU, MC, NU 
PT, SB). 



Published 

With htiernaiionat search report 

Before the expiration of the time limit for amending the 
claims and to be rcpubdshed in the event af the receipt of 
amendments. 



(54) Title: FINITE AMPLrRTDE DISTORTION-BASED INHOMOGENEOUS PULSE ECHO ULTRASONIC IMAGING 
(57) Abstract 

10 

A method and system for Imaging ^ 
a sample. The method includes the steps 
of generating an ultrasonic signal, directing 
llie. signal into a sample, which signal is 
distorted and contains a first order and 
higher order component signals at first 
and higher frequencies respectively. Tlie 
received distorted signal is processed, and 
an image is fonncd, and then displayed, 
from one of the higher order cornponcni 
signals of the received distorted signal. 




FOR TUB PURPOSES OF INFORMATION ONLY 



Codes used to identify Slates party to the PCT on the front pages of pamphlets publishing international applications under the PCT. 



AL 


AlbwU 


ES 


Spain 


LS 


Lesotho 


SI 


Slovenia 


AM 


Antienia 


FI 


Fiulaiid 


LT 


LithuaTiia 


SK 


Slovalcia 


AT 


Austria 




Fruicc 


W 


LuKcmboui^g 


SN 


Senegal 


AU 


AusiraJli 


GA 


Gabon 


LV 


Uivla 


S2 


Swaziland 


A7. 


Azerbaijin 


OB 


United Kb^dom 


MC 


Monaco 


TD 


Chad 


&\ 


Bosnia and lUne^CMm 


GE 


Grmgia 


MD 


Republic of Moldora 


TG 


Togo 


BB 


Bocbadot 


GH 


Ghaiia 


MG 


Marlaga^C2r 


TJ 


TajilusUn 


BE 


Qclgiuoi 


GN 


Guinea 


MK 


'Itc formsr Yugpsliv 


TM 


TUilanenistui 


BI' 


fiurlrina Faio 


GR 


Greece 




Republic ue' Macedonia 


TR 


Turkey 


BG 


Bulgaria 


HU 


Hungaty 


ML 


Mali 


TT 


Tlrinidxd and Totiago 


BJ 


Benin 


IZ 


Treland 


MN 


Mongolia 


DA 


Ukraine 


BR 


BiBZil 


IL 


latast 


MR 


Mauritania 


UG 


Uganda 


BY 


Belarus 


IS 


Iceland 


M\V 


Malawi 


US 


Unhed States of America 


CA 


Canada 


IT 


TUily 


MX 


Mexico 


U7. 


llzbelnstan 


CP 


Cenuni African Republic 


jr 


Japan 


KE 


Nigcx 


VN 


Via Nam 


CC 


Conga 


KK 


Kenya 


NL 


Netherlands 


YU 


Yugoslavia 


ClI 


SwUzerluid 


KG 


Kyrg)'2stan 


NO 


Norway 


zw 


Zimbabwx 


CI 


C&te d'Woirc 


KP 


Democretu: People' a 


St 


New Zealand 






CM 


Can)en>on 




Republic of Korea 


PL 


Poliind 






CN 


China 


KR 


Republic of Korea 


PT 


Portugal 






cu 


CuU 


KZ 


Kozalcston 


RO 


Romania 








Czech Itepublic 


LC 


jfaint Luda 


RU 


Russian t«dcxalion 






OE 


Cennany 


U 


Ltedueuatein 


SD 


Sudan 






DK 


Denmark 


LK 


Sri Litnka 


SE 


Sweden 






ILK 


Kstonia 


LR 


T^ibeiia 


SG 


Siagapcnt 







wo 98/20361 



PCT/US97/20162 



0 



-1- 



PTNTT^P fiMPLITUPy ; nTfiT^n^TTQN- BASED 
INH0MQQENEQU5 PULSE ECHO U LTRASONIC IMAGING 

5 This invention generally relates to 

ultrasonic pulse echo imaging, and more specifically, 
to ultrasonic pulse echo imaging based on the 
distortion of ultrasonic signals transmitted into 
samples. 

10 Ultrasonic pulse echo imaging is widely used 

in many medical applications. While this technique 
has received wide acceptance, it would be desirable to 
improve the resolution of the Images formed from this 
technique. For example, ultrasonic pulse echo imaging 

15 in inhomogeheous media suffers from significant 
lateral and contrast resolution losses due to the 
defocusing effects of the inhomogenei ties of the 
media. The losses in lateral and contrast resolution 
are associated with increases in the width of the main 

20 beams and increases in side lobe levels, respectively. 

These two forms of resolution loss represent 
a significant hurdle Lo improving the clinical utility 
of biomedical ultrasonic imaging. A number of 
research efforts are currently underway to investigate 

25 the defocusing effects of tissue and to consider 
corrective measures. These efforts, however, 
generally assume linear propagation and base the 
image -formation process on the reception of the 
transmitted pulse, 

30 
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The present: invention is directed to 
1 improved ultrasonic imaging and methods, as well as 
improved resolution of images formed from ultrasonic 
echo signals. The present invention utilizes the 
enhanced inhomogeneous focusing properties of the 
5 finite amplitude distortion generated higher harmonics 
of an ultrasonic imaging beam in order to obtain 
improved contrast resolution and lateral resolution 
jinngns. 

More specifically, the present invention is 
directed to a method which includes the steps of 
generating an ultrasonic 'signal , directing the signal 
into a sample, which signal is distorted and contains 
a first order and higher order component signals at 
first and higher frequencies respectively. The 
^5 received distorted signal is processed, and an image 
is formed, and then displayed, from one of the higher 
order component signals of the received distorted 
signal . 

With the preferred embodiment of the 
20 invention disclosed herein in detail, the ultrasonic 
image is based on one of the received finite amplitude 
distortion component (or nonlinearly -generated higher 
harmonics) associated with the transmitted signal. In 
the simplest case, in which the transducer emits 
25 negligible energy in the second harmonic bandwidth, 

such an image can be formed by adding an initial high 
pass filtering of the received signal. In general, 
such an image can be formed by using a two pulse 
transmit, receive, normalize, and then high pass 
30 filtering scheme. Such a two pulse scheme can be used 
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to remove the source or linear content in the second 
^ and higher harmonic bandwidths. 

Further benefits and advantages of the 
invention will become apparent from a consideration of 
the following detailed description given with 
5 reference to the accompanying drawings, which specify 
and show preferred embodiments of the invention. 

Figure 1 shows an ultrasonic imaging system 
embodying the present invention. 

Figures 2a, 2b and 2c show various 
'^^ parameters associated with the linear propagation 
results for a focused 2 MHz Gaussian Transducer 
operating in a liver medium. 

Figure 3 shows discrete harmonic velocities 
used to compute various values associated with a 2 MHz 
propagation in a liver medium. 

Figures 4a and 4b display nonlinear 
propagation results for a focused 2 MHz Gaussian 
transducer. 

Figure 5a shows the log- scaled, normalized 
one-way focal plane profiles of the 2 MHz fundamental, 
the 4 MHz second harmonic, and the 4 MHz fundamental, 
and Figure 5b shows the corresponding two way profiles 
for these beam patterns. 

Figure 5a shows the on-axis source plane and 
25 the subsequent focal plane for a 2 MHz Gaussian 
source. 

Figure 6b illustrates the corresponding 
spectrums of the source and focal planes shown in 
Figure 6a. 

30 Figure 6c shows the focal nonlinear 
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distortion pulse obtained by constructing the waveform 
1 using only the focal spectral information shown in 
Figure 6b from 3 to 8 MHz. 

Figure 7 shows the focal plane profile from 
a 2 MHz continuous v/ave propagation and a 2 MHz 
5 centered pulse propagation. 

Figure 8 is a table giving the on -axis 
source plane intensity versus the received second and 
third harmonic levels and focal field parameters. 

Figure 9 is a Table providing focal and 
received second harmonic levels versus focal length. 

Figures 10a and 10b show the log- scaled 
first, second and third harmonic axial amplitudes for 
the focused 2 MHz Gaussian transducer and the 
corresponding log- scaled focal plane radial beam 
15 profiles , 

Figure 11 illustrates scaled axial 
amplitudes for a 4 MHz second harmonic and 4 MHz 
fundamental harmonic beams in a liver medium. 

Figures 12a- 12f depict pairs of one-way 
20 focal plane harmonic amplitude diameters and the 
corresponding average radii for an unjittered, or 
homogeneous, path and an abdominal wall - j ittered 
propagation path. 

Figures 13a -13d show two-way average radii 
25 results for the abdominal wall- jittered propagation 
path represented in Figures 12a- 12f. 

Figures 14a and 14b display normalized two- 
way averaged radial results from five abdominal wall- 
j ittered propagations for the 2 MHz fundamentals, the 
30 4 MHz second harmonics and the 4 MHz fundamentals, and 
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the corresponding radially integrated magnitudes. 
1 Figures 15a and 15b show normalized two-way 

averaged radial results from five abdominal wall- 
jittered propagations for the 4 MHz fundamentals, the 
8 MHz second harmonics and the 8 MHz fundamentals, and 
5 the corresponding radially integrated normalized 
magnitudes. 

Figures 16a and 16b show normalized two-way 
averaged radial results from five breast jittered 
propagations for the 2 MHz fundamentals, the 4 MHz 
second harmonics and the 4 MHz fundamentals, and the 
corresponding radially integrated magnitudes. 

Figure 17 is a table giving the -20dB full- 
widths for the average two-way profiles shown in 
Figures 14a, 15a and 16a and the full -widths at the 
15 0.9 level for the integrated profiles of Figures 14b, 
15b and 16b. Also shown are the corresponding results 
from the 8 MHz breast jittered propagations. 

Figures 18a -18d show an imperfect source 
pulse, linearly- scaled and log- scaled focal spectruras 
2^ of that pulse, and the corresponding nonlinear 

distortion pulse obtained by constructing with the 
spectral information in Figure 18c. 

Figure 19a shows the log -scaled focal 
spectrum of Figure 19c overlaid with the focal 

25 

spectrum from the same source using a half amplitude 
version of the source pulse, as depicted in Figure 
19c. 

30 Figure 19b shows the resulting difference 



35 



wo 98/20361 PCT/US97/2016Z 

-6- 



spectrum computed for the tv;o spectrum shown in Figure 
1 20a. 

Figure 19c shows the corresponding nonlinear 
distortion pulse obtained by constructing with the 
spectral information in Figure 19b and starting at 
5 2.75 MHz. 

Figure 1 illustrates ultrasonic imaging 
system 10. A pulse generator 12 and a function 
generator 14 produce a sinusoidal pulse ultrasonic 

10 signal of, for example, 2.0 MHz at a pulse repeat 
frequency of, for instance, 1 kHz. This signal is 
sent to amplifier 20, which amplifies the signal and 
transmits the amplified signal to transducer- receiver 
unit 22, and this unit then transmits the signal into 

15 sample 24. 

m this sample 24, the input signal is both 
distorted and reflected. The distortion creates a 
distorted signal having a multitude of component 
signals, each of which has a respective frequency or 

20 frequency bandwidth. The distorted signal is 

reflected by sample 24. and this reflected signal is 
received by transducer -receiver unit 22, amplified by 
pre- amplifier 30, and then further amplified by 
amplifier 32. The received and amplified signal is 

25 then sent through a high-pass filter 34 to enhance the 
relative strength of the desired higher harmonic 
component of the received signal. The resulting 
signal is digitized in analog- to-digital converter 36, 
and then processed by processor 40 Lo produce an 

30 image. 
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This image may be displayed on a video monitor 

1 42, stored on a video cassette recorder (VCR) 44, 

output on a printer device 46, or stored in any of a 
variety of hard copy storage devices 50, such as 
medical film recorders, digital tape machines, optical 
5 disks, magnetic tapes and disks or the like. Suitable 
means may be used to move the focal point of the 
transmitted signal around sample 24. For instance, 
transducer -receiver unit 22 may be a phased array unit 
having electrical circuitry to move the focal point of 
the transmitted signal. Alternatively, a motor 52 may 
be employed to move transducer -receiver unit 22 and 
thereby move the focal point of the transmitted signal 
around sample 24. 

System 10 can also be operated in a two 

^5 pulse scheme or mode. In this mode of operation, 

system 10 generates and transmits into sample 24 two 
different pulses. Preferably, the transmitted signals 
are identical except that one of them is scaled up in 
pressure. The pulses are transmitted one after 

20 another with, for example, approximately 1/4000 second 
interval between them. The reflected, distorted 
signals from both pulses are received by transducer 
unit 22, and these signals are then digitized in 
analog- to-digital converter 36. The digital data 

25 values obtained from the first pulse are stored, 
scaled and then subtracted from the digital data 
values obtained from the second pulse, producing a 
difference or resultant signal. Subsequent processing 
of this difference signal in system 10 is the same as 

30 in the above- described one pulse case. One suitable 
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procedure for scaling the digital data values obtained 
from the first pulse is discussed below. 

Several computational models exist which can 
accurately describe the finite amplitude propagation 
of a continuous beam. Such models can be extended to 
5 compute the case of a propagating pulse as well. 

These models account for the effects of diffraction, 
nonlinearity and absorption. One of these models, 
referred to as the NLP model is described in "^3ew 
approaches to nonlinear diffractive field propagation" 

10 J. Acoust. Soc. Am. 90, 488-499 (1991) by P.T. 

Christopher and K.J. Parker, the disclosure of which 
is herein incorporated by reference, and this model 
was used to compute the following linear and nonlinear 
examples- This model has been updated to account for 

15 the effects of dispersion. The associated harmonic 
velocities were computed using an algorithm disclosed 
in "Modeling acoustic field propagation for medical 
devices," Ph.D. thesis. U. of Rochester (1993) by T. 
Christopher, the disclosure of which is herein 

20 incorporated by reference. 

The propagations discussed below are for a 
Gaussian apodized, axial ly symmetric focused source. 
This form of device offers excellent image quality and 
produces a field comparable to that of the two 

25 dimensional array-based transducers now being 
developed. 

Many biomedical imaging devices are not 
axially symmetric, though. The finite amplitude beams 
produced by such devices are well described by the 
30 computations for the axially symmetric transducers 
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discussed here in detail, though. In measuring the 
nonlinear harmonic generation from an unfocused 
rectangular source, Kamakura, Tani, Kumamoto and Ueda 
noted in "Harmonic generation in finite amplitude 
sound beams from a rectangular aperture source," ON 
AcQush. 5oc. ;^jn, 97. 3510-3517 (1995), that "the 
[nonlinear] harmonic pressure levels in the far field 
[v/ere] almost the same as from a circular aperture 
source with equal face area and equal initial 
pressure, independent of the source levels." Though 
this result was obtained for only one device (with a 
ratio of source side lengths or aspect ratio of 11 to 
6), the higher harmonic pressure levels associated 
with a non-axially symmetric device are approximately 
equal to those of the corresponding axially symmetric 
source. More importantly, the lack; of axial siTiimetry 
does not affect the relative sidelobe advantages 
exhibited by the nonlinear ly- generated harmonics in a 
homogeneous propagation. These homogeneous path 
sidelobe level advantages are the basis for the 
imaging -relevant advantages of the higher harmonics in 
an iiiliomogeneous propagation. 

Figures 2a, 2b and 2c show various 
parameters associated with the linear propagation 
results for a focused 2 MHz Gaussian transducer 
operating in a liver medium. In particular, the 
source plane amplitude profile, the on- axis amplitude, 
and the radial focal plane (at Z = 6cm) beam profile 
are shown in Figures 2a, 2b and 2c respectively. The 
results shown in Figures 2a, 2b and 2c were obtained 
by computing the linear, liver path propagation of the 
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field of a focused 2 MHz Gaussian source using the NLP 
beam propagation model. The NLP model . propagates a 

1 planar, normal velocity description of the acoustic 
field. No inhomogeneiUies or phase aberrations were 
accounted for in this propagation or any of the 
subsequent propagations considered immediately below, 
5 The relevant liver propagation parameters used were c 
= 1570 m/s, p = 1.05 g/cm\ = 0.03 Np/cm and b = 1.3 
- (where « and b are the coefficients describing 
absorption in a power law form) . 

The Gaussian shading of the magnitude of the 

10 source plane normal velocity field was such that the 
half -amplitude radial distance was 0,84 cm. The on- 
axis, source plane RMS acoustic intensity (pcIu|V2, 
where u is the acoustic particle velocity) for the 
field was 2 W/cm^ The radial extent of the source 

15 was 1.5 cm. The source plane field was focused using 
a spherically- focusing factor (e^®'''\ where e(r) = 
(2nf /c)/P+P) - The geometric focal length F was 6 cm 
and the sound speed (c) used to compute e(r) was that 
of water (1500 m/s) . 

20 Figures 2a and 2b depict the normal velocity 

magnitudes of the Gaussian transducer's source plane 
and on -axis fields respectively. Figure 2c displays 
the focal plane (z = 6 cm) radial profile of the 2 MHz 
field. The drop in the magnitude of the field from 

25 the mainlobe to the first sidelobe in Figure 2c is 36 
dB. In the absence of strong medium phase aberration 
this allows the device to produce high contrast 
images , 

The same 2 MHz Gaussian source was then 
30 propagated nonlinearly through the same liver path. 
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The nonlinear parameter |5 used to represent liver was 
4.7. The NLP model used 4 harmonics (2, 4, 6 and 8 

1 MHz) to compute the pre -focal region (z = 0 to z = 3 
cm) propagation and up to 10 harmonics to represent 
the subsequent focal and post-focai region 
propagation. The harmonic velocities were computed 

5 using the algorithm described in "Modeling acoustic 
field propagation for medical devices." T. 
Christopher, Ph.D. Thesis, University of Rochester 
(1993) . The fundamental or 2 MHz component had a 
propagation speed of 0.157 cm/microsecond (given above 
10 as c) . The discrete harmonic velocities used by NLP 
to compute the 2 MHz propagation are shown in Figure 
3. 

Figure 4a displays the axial magnitudes of 
the fundflmental. second harmonic, and third harmonic 

15 fields, at 102a. 102b and 102c respectively, as 
computed for the nonlinear propagation. The 
fundamental or 2 MHz axial curve is only slightly 
different from the corresponding linear curve shown in 
Figure 3b. At z = 6 cm the nonlinear 2 MHz curve is 

20 about 1% lower than the 2 MHz linear curve. This 

difference was due to growth of the higher harmonics 
in the nonlinear propagation. In Figure 4b the 
corresponding focal (z = 6 cm) pressure waveforms from 
the linear and nonlinear computations are displayed at 

25 104a and 104b respectively. The pressure waveforms 
were obtained by converting NLP's normal velocity 
output to pressure using the impedance relation (all 
subsequent pressure waveforms were obtained in this 
way) . 

20 The 2 MHz or fundamental beam pattern 104b 
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associated with the nonlinear propagation is almost 
identical with the 2 MHz beam pattern 104a of the 
1 linear propagation. Only in a linearly scaled overlay 
plot of the two beam patterns are there visible 
differences. These differences are very small and are 
limited to the near axis portions of the beam 
5 patterns. Only at much higher source intensities are 
the effects of nonlinearity significant to the details 
of the fundamental's field. These results are 
consistent with the empirical observation that linear 
modeling of biomedical ultrasonic device fields 
10 accurately describes their (linear -based homogeneous 
path) imaging performance. 

Figure 5a depicts the 2 MHz fundamental and 
4 MHz second harmonic focal plane beam amplitude 
profiles at 106a and 106b. Also shown in Figure 5a at 
15 106c is the corresponding 4 MHz fundamental profile. 
The 4 MHz fundamental result was obtained by computing 
the linear propagation of the same Gaussian transducer 
operating at a source frequency of 4 MHz. All three 
beam profiles in Figure 5a have been normalized to 
20 have on -axis field magnitudes of 1. The finite 

amplitude distortion- generated second harmonic focal 
profile 106b has a slightly broader mainlobe than the 
corresponding (4 MHz) fundamental profile lC6c. The 
radial half -amplitude distance of the second harmonic 
25 profile is 36% greater than that of the 4 MHz 

fundamental profile (0.0983 cm versus 0.0723 cm). The 
second harmonic profile also has much lower sidelobes 
than the 4 MHz fundamental profile. 

For imaging purposes, the two-way focal 
30 plane beam pattern of the Gaussian transducer is of 
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interest. The two-way focal beam pattern accounts for 
both the characteristics of the transmitted pulse in 
1 the focal plane and the corresponding sensitivity of 
the transducer to pulses reflected back from this 
plane. For linear propagations, the two-way beam 
pattern for a given depth can be obtained by squaring 
5 the corresponding transmit or one-way beam pattern. 
In Figure 5b the normalized two-way linear beam 
patterns for the Gaussian transducer operating at 2 
and 4 MHz are depicted at 110a and 110b. These curves 
were obtained by squaring the corresponding one-way or 

10 transmit beam patterns shovm in Figure 5a. 

Also shown in Figure 5b at 110c, is the two- 
way focal plane beam pattern associated with the 4 MHz 
second harmonic field. Since the amplitudes of the 
reflected pulses are much smaller than the transmitted 

15 pulses, the propagation of the reflected field back to 
the transducer is essentially linear* Thus, the two- 
way focal plane beam pattern shown at 110b for the 
second harmonic was obtained by multiplying the 
corresponding one-way pattern 106b shown in Figure 5a 

20 by the 4 Mii2 fundamental one-way pattern 106b also 
shown in Figure 5a. 

The second harmonic's two-way beam pattern 
has a half -amplitude mainlobe width (or -6 dB 
beamwidth) that is 12% greater than that of the 

25 corresponding beam pattern 106c, The -20 dB beamwidth 
of the second harmonic is 13% greater than that of the 
4 MHz fundamental. The sidelobe advantage displayed 
in the focal plane profiles of Figure 5a is maintained 
in the two-way results shown in Figure 5b, 

These homogeneous results show that the 
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second harmonic field of a focused, apodized 
transducer offers advantages in contrast resolution 
1 over the corresponding linear transducer field. 

A pulse propagation was next considered for 
the 2 MHz Gaussian source. The on-axis, source plane 
pressure pulse used is displayed in Figure 6a as the 
5 desired curve 112a. This pulse was computed by 
applying a Gaussian window to a 2 MHz cosine. The 
peak pressure of the pulse was the same as for the 2 
W/cm^ continuous case considered above. The radial 
amplitude shading or apodization and the spherical 

10 focusing ot the source field were also the same as in 
the previous continuous wave case. The initial source 
plane pulse consisted of 128 samples across 8 
microseconds , 

The magnitude of the Fourier transform of 

15 the 8 microsecond long source pulse .is shown at 114a 

in Figure 6b. A straightforward implementation of the 
nonlinear imaging system and method requires 
negligible overlap between the sources 's spectral 
bandwidth and that of the nonlinear second harmonic 

20 {more generally, this would also insure negligible 
overlap between any of the successive harmonic 
spectral bands) . The source spectrum 114b shown in 
Figure 6b meets this requirement. Nonlinear images 
based on source pulses v/ith broader spectrums or with 

25 significantly more energy in the second harmonic 

bandwidth than the one depicted at 112b and 114b in 
Figures 6a and 6b could be obtained by using an 
alternative nonlinear imaging scheme described below. 

The source plane was then defined using the 

30 64 harmonic Fourier transform of the source pulses- 
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This multiharmonic source radius was then input into a 
pulse-propagating version of the NLP model (a model 
1 similar to the lithotripter model presented in 
''Modeling the Dornier HM3 Lithotripter." T. 
Christopher. J. Acoust, 5oc. Am. 3088 • 3095 (1994). 
The focal output of the resulting nonlinear pulse 
5 propagation is also shown in Figures 6a and 6b as 
solid curves. The focal pulse waveform has slightly 
smaller peak positive and negative pressures than the 
corresponding continuous waveform shown in Figure 4b. 
Consistent with the smaller amplitudes, the focal 

10 pulse is also less distorted than the continuous 

waveform. The ratio of the second harmonic's focal 
amplitude to that of the fundamental's for this pulse 
propagation was 70% of the same ratio for the 
corresponding continuous source considered above. 

15 In Figure 6c, the waveform associated with 

only the spectral bands of the second, third, and the 
first half of the fourth harmonic (3 through 8 MHz) is 
shown. This waveform was computed as a high pass 
filtered reconstruction of the spectral information 

20 depicted in Figure 6b, A rectangular window with a 
transition at 3 MHz was used in filtering the 
transform data. The straightforward nonlinear imaging 
approach disclosed herein may use a distortion pulse 
like that shown in Figure 6c in order to image the 

25 scattering medium. 

Not shown for the pulse propagation 
considered are the axial and radial harmonic 
descriptions. These descriptions were found to be 
identical in form to those computed for the nonlinear 

30 propagation of the continuous 2 MHz field. All other 
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Gaussian-windowed cosine pulses were found to produce 
identical axial and radial harmonic patterns- The 
1 temporal peak amplitude profile of the high pass 

filtered- focal plane data also very closely followed 
the form of the 4 imz second harmonic's amplitude 
profile. The temporal peak amplitude profile for the 
5 pulse propagation considered herein is shown at 116a 
in Figure 7, This profile was computed using the 
focal plane frequency data in the range of 3 to 8 MHz 
(like the waveform in Figure 6c) . Also shown at 116b 
and 116c in Figure 7 are the corresponding 4 MHz 

10 second harmonic profile and a second harmonic bandpass 
filtered (3 to 5 MHz) - temporal peak amplitude profile 
for the same pulse propagation. The similarity of the 
two peak amplitude profiles suggests that there is 
very little energy above 5 MHz for this pulse 

15 propagation. The nonlinear harmonic beam and peak 
amplitude patterns produced by any reasonable source 
pulse can be computed using an appropriate continuous 
approximation of the source. 

In order for the nonlinearly -generated 

20 higher harmonics to be available for imaging in an 
inhomogeneous media, the received higher harmonics 
amplitudes cannot be too low relative to the received 
linear (or transmitted) signal and the transducer's 
dynamic range. Also, to be available for imaging use, 

25 preferably the associated in vivo field amplitudes 
have a mechanical index less than 1.9, 

For the previously considered nonlinear 2 
MHz, liver-path propagation, the amplitudes at the 
focus (z = 6 cm) of the fundamental and second 

30 harmonic were 0.943 MPa and 0.166 MPa, respectively. 
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The second harmonic amplitude was thus 15.09 dB below 
the amplit)ide of the 2 MHz transmitted wave. As a 
1 wave reflected at the focal point travels the 6 cm 
back to the transducer/receiver, frequency dependent 
attenuation reduces the second harmonic by an 
additional 5.62 dB relative to the 2 MHz component. 
5 The received 4 MHz component would thus be 20.71 dB 
below the transmitted 2 MHz component. The 
corresponding figures for the 6 MHz third harmonic are 
27 dB down at the focus and 39,23 dB down for received 
signals. This calculation, it should be noted, does 

10 not include the effect of the possible additional 
relative weakening of the second harmonic received 
signal due to coherent reflectors at the focus in 
combination with the smaller mainlobe o£ the second 
harmonic. This effect is not significant to many bio- 

15 ultrasound imaging applications, though. 

Table 1 of Figure 8 displays the second and 
third harmonic received levels at on-axis, source 
plane intensity values of 0,5, 1, 2, 4 and 8 W/cm^ 
For a given source plane intensity, pulsed devices 

20 would have slightly larger received level 

differentials than those shown in Table 1, For the 
pulsed propagation considered above, this additional 
gap would be 0.86 dB for the second harmonic (based on 
second harmonic bandpass reconstruction and a 

25 comparison of the received peak positive pressure 
levels) . 

Current biomedical ultrasonic imaging 
transducers have dynamic ranges of about 100 dB. Even 
with decreased sensitivity above the transmit 
30 frequencies, these devices are capable of creating 
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second harmonic images. This capability has been 
demonstrated by the creation of second harmonic 
1 contrast agent -response images. Alternatively, a 
separate receiver device with appropriate frequency 
response in the desired nonlinear distortion bandwidth 
can be used. 

5 The effect of focal length on the received 

second harmonic levels for this Gaussian transducer 
operating at 2 W/cm^ id shown in Table 2 of Figure 9, 
From a focal length of 4 cm up to a focal length of 12 
cm the received second harmonic levels dropped off by 

10 7,32 db. Also shown in Table 2 are the corresponding 
focal second harmonic levels. These levels remained 
very constant and thus revealed the decreases in 
received levels as almost entirely due to increased 
return trip distances. 

15 The peak positive and negative pressures of 

the in vivo nonlinear waveform shown in Figure 4b were 
1.12 and -0.84 MPa, respectively. The -0,84 peak 
rarefaction pressure corresponds to a mechanical index 
(MI) of 0.59. The highest preferred level for the 

20 mechanical index is 1.9. Shown in Table l-are values 
of the computed minimum focal pressure and associated 
mechanical index for this and four other values of 
source plane intensity. The minimum pressures and 
thus mechanical indices given in Table 1 have been 

25 corrected for the effects of nonlinearity . A linear- 
only computation would result in larger negative 
pressures and MI values, in particular at the highest 
two source intensity levels. 

The numbers shown in Table 1 show that for 

30 in vivo propagations similar to the one considered 
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here, finite amplitude distortion-based images are 
readily obtainable v/ithin the current mechanical index 
1 safety limit. Even at the lowest source intensity 
case considered (0.5 W/cm^) , a largely second 
harmonic-based image can be obtained by simply 
filtering out the transmitted frequency or 
5 frequencies. Additionally the second and third 

harmonic received levels offer some real-time feedback 
on the magnitude of the focal field amplitudes 
themselves. Finally, the results displayed in Table 2 
suggest that second harmonic imaging may be available 

10 at a wide range of focal depths. 

The formation of the higher harmonic 
constituent beams in a propagation finite amplitude 
beam is a continuous process. In the case of the 2 
MHz Gaussian- shaded, focused beam considered above, 

15 the 4 MHz second harmonic, the 6 MHz third harmonic, 
and additional higher harmonics are continuously and 
cumulatively produced by the beam as it propagates 
away from the source. Of interest here is the 
production and focusing of these nonlinear higher 

20 harmonic beams between the source and focal plane. 

The origin of the higher harmonic beams is 
the ongoing nonlinear distortion of the propagating 
waves comprising the (total harmonic) focused beam. 
The physical effects of diffraction and absorption 

25 concurrently act on the higher harmonic beams and thus 
further define their propagation as well as contribute 
to changes in the resulting focused beam. The NLP 
model assumes that the nonlinear or finite amplitude 
distortion acts in a plane wave fashion on the waves 

30 comprising the focused beam. The NT.»P model uses the 
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frequency domain solution to Burgers' equation in an 
incremental Az fashion to account for this plane wave 
1 distortion approximation. 

The frequency domain solution to Burgers' 
equation as used in the NLP model can be written 

5 

■5F ^ 2C^ 



Jr = At =s II 



t 



11= 1,2,...,7V 



10 where / is the fundamental frequency and u^ is the nth 
term in an N term complex Fourier series describing 
the temporal normal velocity waveform at a radial 
point in the radial description of the field. The 
first summation term in the parentheses represents the 

15 accretion of the nth harmonic by nonlinear combination 
of other harmonics that have a sum frequency of n/. 
The second summation term (with conjugation) 
represents the depletion- of the nth harmonic to other 
harmonics with a difference frequency of nj. For the 

20 case of a Fourier representation of a (periodic) pulse 
waveform, this accretion and depletion of harmonics 
results in some interesting phenomenon including the 
production of a distortion bandwidth below the 
fundamental bandwidth. 

25 Of interest are the terms in equation (1) 

contributing to the growth of the second harmonic and, 
to a lesser extent, the third harmonic, when n equals 
2 in equation (1) , the positive 

30 contributions to T^'' come from the lu^u^ term comprising 



35 



wo 98/20361 



PCT/US97/20162 



the first summation. The negative contributions to 
the. second harmonic are represented by the second 
1 summation in equation (1) and for the propagations 
under consideration here can be approximated by the 
second term in that summation^ 2u3Ui. Even this term, 
though, is negligible throughout most of the source to 
5 focal region propagation due to the relatively small 
amplitudes of the third harmonic. Thus, the nonlinear 
production of the second harmonic throughout most of 
the relevant propagation region is simply proportional 
to the square of the amplitude of the fundamental. 

10 When n equals 3 in equation (1), the 

positive contributions to |j* come from the first two 
terms comprising the first summation. These tv;o terms 
sum to 2UjU^, The first term of the corresponding 
negative contributions to the third harmonic is 3u4Uj. 

15 This term and it's successors are negligible for all 
but the focal region of the highest intensity 
propagations considered herein. 

The harmonic sources of the second harmonic 
and third harmonic are thus lu^u^ and Su^Uj, 

20 respectively. At a given point in the field of the 

propagating beam then, the finite amplitude production 
of the second harmonic is proportional to the square 
of the fundamental harmonic's amplitude. The 
production of the second harmonic off the beam axis is 

25 very small since the amplitude 'of the fundamental beam 
there is quite small. The third harmonic is produced 
in proportion to the product of the first and second 
harmonics and thus its nonlinear production is even 
more strongly weighted towards the beam axis. 

QO Also of interest is how the production of 
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the second and third harmonics vary with the z 
coordinate. Neglecting the effects of absorption and 
1 approximating the effects of focusing by assuming 
spherically ■ converging wave propagation, the 
amplitudes of the fundamental harmonic's mainbeam 
increases approximately linearly with distance from 
5 the source to the focus. At z = F/2 the on -axis 
amplitudes of the beam are about twice the 
corresponding source plane amplitudes. Following from 
the same approximations, the fundamental's mainbeam 
width at 2 = F/2 is about half its source plane width. 

10 Thus, at z = F/2 the area of the fundamental's 

mainbeam is about 1/4 the corresponding source plane 
area. Together, this relation and the previous 
amplitude relation suggest that the second harmonic 
beam production rate versus z is constant (z < F) , 

15 with rate losses due to diminishing fundamental 

mainbeam area balanced by the concurrent gains due to 
increased fundamental amplitudes. 

The fact that the third harmonic production 
is proportional to the product of the fundamental and 

20 second harmonic amplitude, though, implies that the 
production of the third harmonic is strongly weighted 
towards the focal region. In Figure 10a the log- 
scaled, axial amplitudes of the fundamental, second, 
and third harmonics are displayed at 120a, 120b and 

25 120c respectively. The source was the same focused 2 
MHz Gaussian source considered above. The medium's 
parameters were again those of liver. Consistent with 
the above discussion, the amplitudes of the second 
harmonic exhibit a relatively large gain in its growth 

30 from low near field values to significant focal 
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amplitudes. The third harmonic exhibits an even 
higher gain, approximately duplicating the growth in 
1 gain from the fundamental to the second harmonic. 
Both harmonics, though, display post -focal region 
amplitude declines which parallel those of the 
fundamental . 

5 In Figure 10b the corresponding log- scaled, 

focal plane (z = 6 cm, one way) radial beam profiles 
are displayed at 122a, 122b and 122c respectively. 
The fundamental profile drops 49.8 dB over the 1 cm 
radial range displayed. The second harmonic 

10 approximately squares this dp.cline in dropping 85.3 
dB. The third harmonic then continues the 
relationship in dropping 120.5 dB. These declines 
reflect the second and third harmonic, finite 
amplitude production rates discussed above. At z = 8 

15 cm this relationship between the harmonic beam 
profiles continued to hold. 

In Figure 11, the axial amplitudes of the 4 
MHz second harmonic beam, shown at 124a, are overlaid 
with the corresponding 4 MHz fundamental harmonic 

20 beam, shown at 124b. The on -axis source of both of 
the respective propagations was 2 w/cm^. The axial 
curves have been normalized to be unity at z = 6 cm 
and log- scaled. The 4 MHz second harmonic focal 
amplitude was originally (pre -normalization) 15.0 db 

25 before the corresponding fundamental value. The two 
curves are close through the focal region and then 
depart shortly after z = 6 cm as the 4 MHz linear 
curve rapidly declines. 

The results displayed in Figures 10a, 10b, 

30 and 11 show that the 4 MHz second harmonic beam may be 
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less susceptible to the defocusing effects of near 
field phase aberrations than a 4 MHz fundamental beam. 
1 Since only a fraction of the second harmonic beam 
forms in the near field, only this fraction could be 
redirected or defocused by near field jitter. The 
corresponding 2 MHz fundamental beam, though, would 
5 pass in its entirety through the aberration and suffer 
the consequent defocusing effects including increased 
focal plane sidelobe levels. Secondarily, though, 
these higher 2 MHz sidelobe levels could, in turn, 
increase the off-axis nonlinear production of the 4 

10 MHz second harmonic. 

In order to investigate the effects of 
tissue -based phase aberration on the characteristics 
of linear and nonlinear beams, planes of phase delay 
or jitter were introduced into linear and nonlinear 

15 propagations of the focused Gaussian transducer. 

These phase delay planes were computed using measured 
time delays from 5 abdominal wall layers and 5 breast 
wall layers. The 5 abdominal wall specimens had layer 
thicknesses of 2.5, 2.0, 1.5, 1.5, and 1.0-3.0 cm (a 

20 non-uniform layer) giving an average thickness of 1.9 
cm. The 5 breast wall specimens had layer thicknesses 
of 1.5-2.5, 3.0-3.5, 3,5, 4.0 and 2.0-2.5 cm, giving 
an average thickness of 3.0 cm. The average 
thicknesses of the non- uniform layers were used to 

25 compute the 5 -layer averages. All of the measured- 
abdominal v;all and breast wall layers contained an 
outer skin layer. 

Each of the measured abdominal wall time 
delay planes was converted to an equivalent 2 MHz 

30 phase delay plane. Each of these delay planes was 
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then scaled by 0.5 (i.e., each phase delay was reduced 
by a factor of 2) and then applied twice to a given 
1 beam propagation to represent the cumulative aspect of 
the actual tissue delays. In applying a single delay 
plane, the 2 MHz phase delay values were scaled for 
appropriate application to each harmonic present in 
5 the computed field. The first delay plane was 

encountered by the propagating field at z = 0,5 cm and 
the second plane at 2 = 1.5 cm. Further, subdivision 
c and subsequent ap^plications of the abdominal wall 
delay data did not appear to be necessary since it did 

10 not significantly change the resulting focal plane 
fields. Thus, the 2 delay plane application scheme 
satisfied the thin lens approximation. The breast 
delay planes were likewise applied in two steps, the 
first at z = 1 on and the second z = 2 cm. For both 

15 tissue types, the z placement was selected to 

represent the average slice thickness and also to be 
convenient for the Az step size utilized by the linear 
propagation. 

Figures 12a -12f depict pairs of one-way 

20 focal plane harmonic amplitude diameters and the 
corresponding average radii for an unjittered (or 
homogeneous) path and an abdominal wall- jittered 
propagation path. The propagation parameters of both 
mediums were again those of liver. In Figure I2a the 

25 corresponding focal plane diameters for the unjittered 
and jittered 2 MHz fields are overlaid at 132a and 
132b respectively. In Figure 12b the corresponding 
average radii are shown at 134a and 134b respectively 
for the 2 MHz fields. The average radii were obtained 

30 by averaging the focal plane grid of amplitudes around 



35 



wo 98/20361 



PCTAJS97/2ei62 



-26- 

the axis. The corresponding results for the second 
harmonic 4 MHz are shown at 136a, 136b, 140a and 140b 
1 in Figures 12c and 12d, and for the fundamental 4 MHz 
fields are shovm in at 142a, 142b, 144a, and 144b in 
Figures 12e and 12f . In all of these Figures, the 
unjittered field is displayed using solid curves. 
5 In Figures 12a- 12f, two basic effects of 

jitter are visible. The first is the increased 
sidelobe levels associated with the defocusing of the 
abdominal wall phase delays. This sidelobe effect is 
more prominent for both of the 4 MHz fields than for 

10 the 2 MHz fundamental field. The shorter wavelengths 
of these 4 MHz fields allows for greater de- focusing 
by the phase screen. The second harmonic 4 MHz field 
also gets a sidelobe level increase from the 
corresponding increase in the 2 MHz field's sidelobe 

15 levels. 

The second impact revealed in Figures 12a- 
12 f is that there are small changes in the mainlobes 
of all three harmonic fields. In particular, in the 
diameter Figures 12a, 12c and 12e, a shift in the peak 

20 or center of the jittered lobes can be seen. In the 
corresponding average Figures 12b, 12d and 12f though, 
the impact of the jitter is negligible down to 
approximately 20 dB below the peak on -axis value. 
Thus, the impact of this abdominal wall • j ittering does 

25 not appear to involve any significant broadening of 
the mainlobes. Finally, both of the jittered 4 MHz 
mainlobes show decreases in peak amplitudes which are 
consistent with the increased energy present in the 
sidelobe regions. 

30 Figures 13a -13d show two-way average radii 
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results for the same abdominal wall- jittered 
propagation. Each of the curves shovm was obtained by 
1 radially averaging the corresponding two-way planar 
data and then scaling the on- axis value to unity. The 
2 MHz fundamental, 4 MHz second harmonic, and 6 MHz 
third harmonic average radii are shown at 146a, 146b 
5 and 146c in Figure 13a. The 2, 4 and 6 Mllz 

fundamentals are shown at 146a, 150a and 150b in 
Figure 13b. The average two-way sidelobe levels of 
the second and third harmonic can be seen to be 
significantly lower than those of any of the 

10 fundamental harmonics. In Figure 13c, the two 4 MHz 
average two-way profiles are shown at 152a and 152b. 
The 4 MHz fundamental curve has a slightly narrov/er 
mainlobe (9.5% at 20 dB down from peak) and 
significantly higher sidelobe levels than the 4 MHz 

15 second harmonic profile. Likewise, Figure 13d shows 
the 6 MHz fundamental and 6 MHz third harmonic two-way 
profiles at 154a and 154b. The 6 MHz fundamental 
curve is 10,9% narrower at the -20 dB level and also 
shows higher sidelobe levels than the third harmonic 

20 curve. 

The results shown in Figures 13a- 13d show 
that the second and third harmonics maintain lower 
sidelobe levels than the corresponding fundamental 
harmonics in propagating through abdominal wall.. In 

25 Figures 14a and 14b, radial results obtained by 
averaging across 5 abdominal wall - j ittered 
propagations are shown. In Figure 14a the normalized 
radial averages from the 5 two-way planar amplitude 
data sets are depicted. The average 2 MHz result is 

30 overlaid with the corresponding 4 MHz second harmonic 
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and 4 MHz fundamental curves 156a, 156b and 156c 
respectively. The -20 dB width of the 4 MHz 
1 fundamental profile 156c is 6.6% narrower than the 4 
MHz second harmonic profile 156b. The -20 dB width 
for this and multiple propagatipn -averaged results 
are given in Table 3 of Figure 15. 
5 Each of the two-way average profiles shown 

in Figure 14a was then radially summed. Each of the 
average profiles consisted of 499 radial position- 
magnitude value pairs, {ri,m;) , i = 1,...,499. The 499 
pairs discretely described the averaged of 5 abdominal 
10 wall- jittered two-way profiles over a radial extent of 
1.2 cm. The discrete radial summation of the average 
profiles was defined by 499 radial position -summation 
value (r^,Sj) pairs, where each s^ was defined by 

The first term in the summation involves the (r^.m^) 

20 on -axis magnitude value. 

The resulting radial summation or 
integration profiles are shown in Figure 14b, In 
Figure 14b, the integration profiles for the 2 MHz 
fundamentals, the 4 MHz second harmonics and the 4 MHz 

25 fundamentals are shown at 160a, 160b and 160c 
respectively. Each of these integrated two-way 
profiles was scaled such that the value at an off- axis 
radial distance of 1.2 cm was unity. The elevated 
sidelobes of the 2 and 4 MHz fundamental profiles 

30 cause their summation profiles to rise significantly 
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beyond the radial extent of the mainlobe. This 
additional rise represents the potential for 

1 scattering from the sidelobes to significantly reduce 
the contrast resolution of the image. In Table 3 the 
radial extent at which Lhese integration profiles 
reach the 0.9 level is given. This radial extent is a 
5 measure of the sideiobe's potential to reduce the 
contrast resolution o£ an image. In this case, the 
second harmonic's radial extent is 38% less than that 
of it's 2 MHz fundamental and 63% less than the 
corresponding 4 MHz fundamental radial extent. 

10 Alternatively, the summation profiles 

depicted in Figure 14b offer the percent of the two- 
way field's amplitude inside or outside a given 
radius. For example, 91.7% of the 2 MHz fundamental's 
amplitude, 96.7% of the 4 MHz second harmonic's 

15 amplitude, and 88.6% of the 4 MHz fundamental's 
amplitude fall inside a radius of 0.25 cm. The 
corresponding percentages falling outside of 0.25 cm 
are 8.3%, 3.3% and 11.4%, respectively* Ratios of 
these outside percentages could be useful for 

20 inferring the relative contrasts offered in imaging a 
low scattering or void region of a given size. For a 
void region approximately 0 , 5 cm across, the 4 MHz 
second harmonic of this device might then provide 2.5 
times (8 dB) higher contrast than the 2 MHz 

25 fundamental and 3,5 times (11 dB) higher contrast than 
the corresponding 4 MHz fundamental. 

The previous two sets (at 2 and 4 MHz) of 5 
abdominal wall jittered propagations were repeated at 
twice the source " frequencies . Figures 15a and 15b 

30 depict the corresponding averaged results from these 4 
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MHz nonlinear and 8 MHz linear propagations, in 
Figure 15a Lhe radial average amplitude curves 
1 obtained by averaging the 5 two-way data sets are 
shown. Figure 15a shows that the 8 MHz fundamental 
mainbeam 162c is broader than the corresponding 4 MHz 
fundamental mainbeam 162a. The jitter- imposed lateral 
5 resolution limits have been encountered and in fact 
exceeded at this point. Also note that the 8 MHz 
second harmonic mainbeam 162b is narrower than either 
of the fundamentals 162a and 162c. Figure 15b depicts 
the additional sidelobe corruption of the fundamental 

10 beams and the corresponding increase in the second 
harmonic's relative contrast resolution potential. 

Figures 16a, and 16b shov/ average results 
from 5 propagations through breast wall delay data. 
The normalized two-way average radii for the 2 MHz 

15 fundamentals, the 4 MHz second harmonics and the 4 MHz 
fundamentals are shown at 166a, 166b and 166c 
respectively in Figure 16a; and the corresponding 
radially integrated magnitudes are shown at 170a, 170b 
and 170c respectively in Figure 16b. 

20 In Figures 16a and 16b the results 

considered were from 2 MHz nonlinear and 4 MHz linear 
propagations. In Figure 16a, the 4 MHz mainbeams 166c 
are significantly broadened over the corresponding 
abdominal wall- jittered mainbeams in Figure 14a. The 

25 sidelobe levels in Figure 16a are also higher than 
those in Figure 14a. These Figures show that the 
breast wall layers produced significantly more 
distortion than the abdominal wall layers. The second 
harmonic 4 MHz profile 166b is 8.4% narrower at the 

30 -20 dB level than the fundamental 4 MHz 166c. In 
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Figure 16b the integrated profile of the second 
harmonic 170b has a radial extent at the 0.9 level 
1 which is 48% less than the 2 MHz fundamental's 170a 
and 70% less than the 4 MHz fundamental's 17 0c. Thus, 
in the more distorting breast wall - j ittered 
propagations, the relative adx^antages of the 4 MHz 
5 second harmonic were larger than in the abdominal wall 
propagations and included a slightly narrower 
mainbeam. The 4 and 8 MHz results for breast wall- 
j ittered propagations follow closely the developments 
seen in the 4 and 8 MHz curves of Figures 16a and 16b. 

10 Results from these propagations are included in Table 
3 of Figure 17. 

Finally, in all of the jittered propagations 
considered, the second harmonic mainbeam was narrower 
than the fundamental mainbeam. The limits on the 

15 lateral resolution of the linear harmonics eventually 
put limits on the second and other higher harmonics, 
though • 

The above discussions show the liver -path 
beam patterns for a focused Gaussian- apodi zed 

20 transducer operating at 2, 4 and 8 MHz. The non-phase 
aberrated propagations show that the second and higher 
harmonics formed through finite amplitude distortion 
have much lower sidelobe levels than their fundamental 
harmonic or the corresponding linear fundamentals, 

25 The finite amplitude production of these higher 

harmonic beams allow this sidelobe relationship to 
hold for any focused or unfocused transducer. Pulse 
propagation analysis shows that the higher harmonics 
formed in a propagating pulse-beam can be very well 

30 described by considering the harmonics produced in the 



35 



wo 98/20361 



PCT/US5)7/20162 



corresponding continuous wave propagation. Modeling 
results also show that second harmonic levels 
1 sufficient for imaging purposes can be easily obtained 
within the field amplitude limits of the mechanical 
index. 

The introduction of phase jitter as computed 
5 from measured propagation delays from slices of 
abdominal wall and breast wall causes the sidelobe 
levels of the second harmonic and fundamental beams to 
rise significantly. In all of the aberrated 
propagations considered, the two-way profile of the 

10 second harmonic offered narrower -20 dB mainlobe 

widths and lower sidelobe levels than the fundamental 
beam which produced it. These same second harmonic 
profiles had slightly broader mainlobes at 4 MHz in 
abdominal wall -jitter propagations than the 4 MHz 

15 fundamental profiles. In all other jittered 

propagations considered, though, the second harmonic 
offered slightly narrower mainlobes than the eame- 
frequency fundamental and substantially lower sidelobe 
levels. Thus, second harmonic -based ultrasonic images 

20 offer significant improvement in the lateral component 
of contrast resolution. 

One obstacle to obtaining such images are 
artifacts from source contributions to the higher 
harmonic bandwidths. Figure 18a depicts such an 

25 imperfect source wave. This on -axis waveform is 

depicted in particle velocity units and corresponds to 
a peak pressure of one half of that of the pulse shown 
in Figure 6a. The pulse in Figure 18a has the same 
form as that in Figure 6a with the exception of the 

30 initial zero portion. In Figure 18b, the 
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corresponding source pulse spectrum is shown. 
Significant energy content outside the 2 MHz bandwidth 
1 is visible. The resulting computed focal waveform's 
spectrum is depicted in Figure 18c and can be compared 
to the earlier focal spectrum in Figure 6b. The log- 
scaled depiction in Figure 18g shows that a .simple 
5 high pass filtering of the received spectrum contains 
significant contributions from the source. 

In Figure 18d, the resulting nonlinear 
distortion pulse obtained by high pass filtering the 
spectrum of Figure 18c is shown. This pulse can be 
10 compared to the corresponding distortion pulse shown 
in Figure 6c. In both cases, the cut-off frequency 
was 3 MHZ. In Figure 18d the full 8 microsecond 
period of the computed pulse is shown so that the 
ringing associated with the linear or source content 
15 within the second harmonic bandwidth can be seen. 
This ringing could adversely affect the axial 
component of the contrast resolution and the lateral 
gains associated with the lower sidelobe levels. In 
this case, improvements can be obtained by increasing 
20 the source amplitude so as to boost the second 
harmonic bandwidth levels. 

A two pulse scheme may be used to alleviate 
or to eliminate such problems. In this method, two 
source pulses are sent in place of a single pulse in 
25 the image formation cycle. The two pulses are 

identical in form but one is significantly lower in 
amplitude. The received echo from this lower 
amplitude pulse is then used to remove the linear 
content from the high amplitude pulse. This is 
30 accomplished by subtracting an appropriately- scaled 
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version of the received low amplitude signal from the 
corresponding high amplitude signal. The resulting 

1 difference signal may then be high-pass filtered 
followed by the norma] sequence of image formation 
steps. The high-pass filtering is preferred since 
pulse analysis has revealed that the low frequency 

5 content of the difference signal is radially wide- 
spread and reduces the higher harmonic sidelobe 
advantage. 

Figures 19a. 19b and 19c show results 
relevant to an implementation of the tv/o pulse scheme 

10 using the pulse shown in Figure 19a as the high 
amplitude source pulse. In Figure 19a. the focal 
spectrum from Figure 18c is shown at 176a and is 
overlaid with the focal spectrum 176b produced by a 
half -amplitude version of the same source. The low or 

15 half amplitude spectrum was multiplied by two and 

subtracted from the high amplitude spectrum to obtain 
the difference spectrum depicted in Figure 19b. Note 
that if the propagations had not involved finite 
amplitude distortion, then this difference spectrum 

20 would have been all zero. The difference spectrum was 
then high pass filtered and inverse transformed to 
obtain the effective on-axis distortion imaging pulse 
shown in Figure 19c. The transition for the high pass 
filtering was at 2.75 MHz, This same frequency was 

25 also appropriate for filtering off-axis difference 
spectra. 

This two pulse scheme appears to be capable 
of extracting the desired largely- second harmonic 
images from any realistic ultrasonic imaging pulse, 
.30 The third harmonic bandwidth depicted in Figure 19b, 
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though, does not appear to be separable from the 
second harmonic bandwidth. The nodal depth between 

1 these harmonics is not deep enough. Thus, third 
harmonic (or largely -third harmonic) images do not 
appear to be easily obtainable with this two pulse 
scheme. Also, a two pulse scheme with a 1/4000 second 

5 interval between the respective high and low amplitude 
source pulses, does not have significant artifacts due 
to tissue or transducer motion. 

The use of the second harmonic (plus a small 
contribution from higher harmonics) to form an image 

10 is an independent alternative to the phase correction- 
based schemes which have been and are being examined 
by other investigators as a means for improving the 
contrast performance of biomedical ultrasonic imaging. 

The higher harmonics also offer additional 

15 opportunities for correcting for beam distortion. The 
amplitude of the third harmonics in the focal region 
is strongly affected by the fundamental's amplitude. 
As discussed above, the nonlinear production of the 
third harmonic is proportional to the product of the 

20 amplitudes of the fundamental and second harmonics. 
Thi.q means that much of the third harmonic's 
production occurs in the focal region where beam 
distortion can reduce the amplitude of the fundamental 
and second harmonic. Thus, an iterative scheme may be 

25 used to correct for beam defocusing using the 

amplitude of the received third harmonic for feedback. 

Second harmonic images also provide for a 
means of reducing speckle. In particular, an image 
formed as the sura of a second harmonic image and the 

30 corresponding fundamental image would have less 
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speckle than either of the constituent images. Since 
the second harmonic is twice the frequency of the 

1 fundamental and has a largely constant phase 

relationship with the fundamental in the mainlobe, the 
second harmonics image's speckle pattern would be 
conveniently out of phase with that of the 

5 fundamental. 

While it is apparent that the invention 
herein disclosed is well calculated to fulfill the 
objects previously stated, it will be appreciated that 
numerous modifications and embodiments may be devised 
10 by those skilled in the art, and it is intended that 
the appended claims cover all such. modifications and 
embodiments as fall within the true spirit and scope 
of the present invention. 

15 



25 



30 
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CLAIMS 

1. A method of imaging a sample comprising the 
1 steps of: 

generating an ultrasonic signal? 

directing the ultrasonic signal into a sample, 
wherein the sample reflects the signal; 
5 receiving the signal reflected by said cample, 

which received signal is distorted and contains a 
first order and higher order component signals at 
first and higher frequencies respectively; 

forming an image from one of said higher order 
10 component signals of the received distorted signal, 
including the step of removing from the received 
distorted signal the first order component thereof; 
and 

displaying said formed image, 

15 

2. A method according to Claim 1, wherein the 
removing step includes the step of high-pass filtering 
the received, reflected distorted signal to remove 
therefrom the first order component thereof. 

20 

3. A method according to Claim 1, wherein: 
the generating signal includes the steps of 

generating first and second ultrasonic signals; 
the directing .step includes the steps of 
25 directing the first and second ultrasonic signals into 
the saJTiple; 

the receiving step includes the step of receiving 
any first and second signals reflected and distorted 
by said sample; 
30 the forming step includes the steps of 



35 



wo 98/20361 PCT/US97/20162 



-38- 

i) subtracting the received second distorted 
signal from the rp.ceived first distorted signal to 

1 produce a resultant signal, and 

ii) forming the image from said resultant signal. 

4. k method according to Claim 3, wherein the 

5 first and second signals are identical except that one 
is scaled up in magnitude by a factor x (greater than 
1) relative to the other and the second signal is 
transmitted after the reception of the distorted first 
signal. 

10 

5. A method according to Claim 4, wherein the 
sample (i) distorts the first ultrasonic signal to 
produce a first distorted signal, (ii) reflects the 
first distorted signal, (iii) distorts the second 

15 ultrasonic signal to produce a second distorted 

signal, and (iv) reflects the second distorted signal, 

6. A method according to Claim 5, wherein the 
forming step includes the step of: 

20 scaling the smaller received distorted signal 

(corresponding to the unsealed transmitted signal) by 
the previously used scale factor x; 

next subtracting this scaled signal to produce a 
difference signal essentially without frequency 
25 content in the original transmitted bandwidth; and 

forming the image from one of said higher order 
component signals of the difference signal. 

7. A method of imaging a sample, comprising the 
30 steps of: 
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generating an ultrasonic signal: 
directing the ultrasonic signal into a sample, 
1 wherein the sample reflects the signals- 
receiving the signal reflected by said sample, 

which received signal is distorted and contains a 

first order and higher order component signals at 
5 first and higher frequencies respectively; 

forming an image from one of said higher order 

component signals of the received distorted signal; 

and 

displaying said formed image; 
10 wherein said higher order component signals 

include a second order component, and the forming step 
includes the step of forming the image from said 
second order component. 

15 a. A method of imaging a sample-, comprising the 

steps of: 

generating an ultrasonic signal: 
directing the ultrasonic signal into a sample, 
v/herein the sample reflects the signal; 
20 receiving the signal reflected by said sample, 

which received signal is distorted and contains a 
first order and higher order component signals at 
first and higher frequencies respectively; 

forming an image from one of said higher order 
25 component signals of the received distorted signal; 
and 

displaying said formed image; 

wherein the directing step includes the step of 
maintaining the sample substantially free of any 
30 contrast agent not naturally present in the sample. 



35 



9. A system for imaging a sample, comprising: 
means for generating an ultrasonic signals- 
means for directing the ultrasonic signal into a 

sample, wherein the sample reflects the signal; 

means for receiving the signal reflected by said 
sample, which received signal is distorted and 
contains a first order and higher order component 
signals at first and higher frequencies respectively; 

means for forming an image from one of said 
higher order component signals of the received 
distorted signal, said means for forming the image 
including* means for removing from the received 
distorted signal the first order component thereof; 
and 

means for displaying said formed image; 

10. A system according to Claim 9, wherein the 
means for removing the first order component from the 
received distorted signal includes a high-pass filter 
to filter the received, reflected distorted signal to 
remove therefrom the first order component thereof. 

11. A system according to Claim 9, wherein: 
the means for generating the ultrasonic, signal 

includes means for generating first and second 
ultrasonic signals; 

the means for directing the ultrasonic signal 
into the sample includes means for directing the first 
and second ultrasonic signals into the sample? 

the receiving means includes means for receiving 
any first and second signals reflected and distorted 
by said sample; 



wo 98/20361 



PCT/US97n0162 



the means for forminq the image includes 

i) means for subtracting the received second 
1 distorted signal from the received first distorted 

signal to produce a resultant signal, and 

ii) means for forming the image from said 
resultant signal, 

5 

12. A system according to Claim 11, wherein the 
first and second signals are identical except that one 
is scaled up in magnitude by a factor x (greater than 
1) relative to the other and the second signal is 

10 transmitted after the reception of the distorted first 
signal . 

13. A system according to Claim 12, wherein the 
sample (i) distorts the first ultrasonic signal to 

15 produce a first distorted signal, (ii) reflects the 
first distorted signal, (iii) distorts the second 
ultrasonic signal to produce a second distorted 
signal, and (iv) reflects the second distorted signal. 

20 14. A system according to Claim 12, wherein the 

forming means includes: 

means for scaling the smaller received distorted 
signal (corresponding to the unsealed transmitted 
signal) by the previously used scale factor x; and 

25 next for subtracting this scaled signal to 

produce a difference signal essentially without 
frequency content in the original transmitted 
bandwidth; and 

means for forming the image from one of said 

30 higher order component signals of the difference 
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signal . 

1 15. A system for imaging a sample, comprising: 

means for generating an ultrasonic signal; 
means for directing the ultrasonic signal into a 
sample, wherein the sample reflects the signal? 
5 means for receiving the signal reflected by said 

sample, which received signal is distorted and 
contains a first order and higher order component 
signals at first and higher frequencies respectively; 
means for forming an image from one of said 
10 higher order component signals of the received 

distorted signal; wherein said higher order component 
signals include a second order component, and the 
means for forming the image includes means for forming 
the image from said second order component; and 
15 means for displaying said formed image. 

16, A system for imaging a sample, comprising: 
means for generating an ultrasonic signal; 
means for directing the ultrasonic signal into a 
20 sample, wherein the sample reflects the signal; 

means for receiving the signal reflected by said 
sample, which received signal is distorted and 
contains a first order and higher order component 
signals at first and higher frequencies respectively; 
25 means for forming an image from one of said 

higher order component signals of the received 
distorted signal; and 

means for displaying said formed image; 
wherein the sample is substantially free of any 
30 contrast agent not naturally present in the sample. 



35 



wo 98/20361 



PCT/US97/20162 



-A3- 

17. A method of imaging a sample, comprising the 
steps of: 

1 generating a series of ultrasonic pulse signals; 

directing the ultrasonic pulse signals into a 
sample, wherein the sample reflects the pulse signals; 
receiving the pulse signals reflected by said 
5 sample, which received pulse signals are distorted and 
contain a first order and higher order component 
signals at first and higher frequencies respectively; 

forming an image from one of said higher order 
component signals of the received distorted pulse 
10 signals; and 

displaying said formed image. 

18, A method of imaging a biological sample, 
comprising the steps of: 

generating an ultrasonic signal; 
directing the ultrasonic signal into the 
biological sample, wherein the sample reflects the 
signal; 

receiving the signal reflected by said sample, 
which received signal is distorted and contains a 
first order and higher order component signals at 
first and higher frequencies respectively; 

forming an image from one of said higher order 
component signals of the received distorted signal; 
and 

displaying said formed image. 



15 



20 



25 



19. A system for imaging a sample, comprising: 
means for generating a series of ultrasonic pulse 
signals; 

30 

means for directing the ultrasonic pulsG signals 
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into a sample, wherein the sample reflects the pulse 
signals; ' 

means for receiving the pulse signals reflected 
by said sample, v/hich received pulse signals are 
distorted and contain a first order and higher order 
component signals at first and higher frequencies 
respectively; 

means for forming an image from one of said 
higher order component signals of the received 
distorted pulse signals; and 

means for displaying said formed image. 

20. A system for imaging a sample, comprising: 
means for generating an ultrasonic signal; 
means for directing the ultrasonic signal into a 

sample, wherein the sample reflects the signal; 

means for receiving the signal reflected by said 
sample, which received signal is distorted and 
contains a first order and higher order component 
signals at first and higher frequencies respectively; 

means for forming an image from one of said 
higher order component signals of the received 
distorted signal; and 

means for displaying said formed image; 
wherein the sample is a biological sample. 

21. A method of imaging a sample, comprising the 
steps of: 

generating an ultrasonic signal; 

directing the ultrasonic signal into the 
biological sample, wherein the sample linearly 
reflects the signal; 
3^ receiving the signal linearly reflected by said 
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sample, which received signal is distorted and 
conLains a first order and higher order component 
1 signals at first and higher frequencies respectively; 
forming an image from one of said higher order 
component signals of the received distorted signal; 
and 

5 displaying said formed image. 

22. A system for imaging a sample, comprising: 
means for generating an ultrasonic signal; 
means for directing the ultrasonic signal into a 

10 

sample, wherein the sample linearly reflects the 
signal; 

means for receiving the signal linearly reflected 
by said sample, which received signal is distorted and 
contains a first order and higher order component 
signals at first and higher frequencies respectively; 

means for forming an image from one of said 
20 higher order component signals of the received 
linearly reflected, distorted signal; and 

means for displaying said formed image. 
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